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Abstract 

We present a linear spin and orbital wave theory to account for the spin 
and orbital orderings observed experimentally in undoped manganite. It is 
found that the anisotropy of the magnetic structure is closely related to the 
orbital ordering, and the Jahn- Teller effect stabilizes the orbital ordering. 
The phase diagram and the low energy excitation spectra for both spin and 
orbital orderings are obtained. The calculated critical temperatures can be 
quantitatively comparable to the experimental data. 
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LaMnOa is the parent compound of colossal magnetoresistance (CMR) manganites, and 
has been studied both experimentally and theoretically. The compound is an insulator with 
layered antiferromagnetic (A-type AF) spin ordering and an orbital ordering of e g electrons.^ 
Murakami et a/.i has recently succeeded in detecting the orbital ordering in LaMnOs by 
using resonant x-ray scattering techniques with the incident photon energy tuned near the 
Mn K-absorption edge. The orbital order parameter decreases above the Neel temperature 
Ttv ~ 1A0K and persists until To ~ 780K. Theoretically, the problem of orbital degeneracy 
in a d-electron system was pioneered by Kugel and Khomskiii in 1970's and investigated 
extensively in recent years. 

In this paper, starting from an effective Hamiltonian of the spin and orbital interactions, 
as well as the JT coupling between the e g electrons and the lattice distortion, we investigate 
the interplay among the spin, orbit and the lattice distortion. We present the phase diagram 
as functions of interaction parameters, and obtain the low-energy excitations of the system 
in different phases. It is found that special properties of the orbital operators can result 
in an anisotropy of the magnetic structure and an energy gap of the orbital excitations. 
We also estimate the critical temperatures for spin and orbital orderings as well as their 
dependence on the JT coupling. The calculated results are comparable to the experimental 
measurements. 

The effective spin and orbital interactions are derived by the projection perturbation 
method up to the second ordeJHl 



ij ij 

-J 3 BSrS i + 6Kn«, (1) 

ij 

where S, is the spin operator of S = 2. The three terms describe three processes 
with different intermediate states. §0 n f = d\ a d ia and nf = 4 & d m are the particle 
number operators of e g electron in orbit states \a) = cos(<p a /2)\z) + sin((p a /2)\z) and 
| a) = — sm(ip a /2)\z) + cos(y? a /2)|z), respectively, with orbital states \z) oc (3z 2 — r 2 )/v^3 
and \z) oc x 2 — y 2 . Here (p a depends on the direction of the (ij) bond by (p x = — 27r/3, 



(p y = 2tt/3, and (p z = 0, respectively, for bond (ij) parallel to the x, y and z directions. 
The introduced d\ a ,d ia and d\ a ,d ia are operators in the orbital space, with d\ a \0) = \a), 
d\&\fy = 1^)5 they should satisfy the constraint n" + n™ = 1. 
The JT interaction can be expressed as0 

Hjt = -9 E 4 7 T 7 y ■ QAy + y E Qi ' ( 2 ) 

177' i 

where T = (T X ,T Z ) are the Pauli matrices in the orbital space with 7 (7') = z or z, 
and g is the coupling between the e g electrons and the local JT lattice distortion Qj = 
Qi(sin0j, cos0,). Here we have neglected the terms for the anharmonic oscillation of JT 
distortion and the higher order coupling, their effect being regarded approximately as giving 
a preferable direction fa of the JT distortion observed experimentally 

Experimental measurement on LaMnOs indicates that the critical temperature of the 
orbital ordering, To, is much higher than that of the magnetic ordering, T/v,i As a result, 
the spin and orbital degrees of freedom, which are coupled to each other in Hamiltonian (0), 
may be separately treated by the Hartree-Fock mean-field approach. The total Hamiltonian 
is reduced to H MF = H s + H + E , where H s and H Q are the decoupled spin and orbital 
Hamiltonians, respectively, and E Q is an energy constant. The spin Hamiltonian Hs is given 
by 

Hs = 2J Jifii ' Sj , (3) 

m 

with the effective spin coupling depending on the orbital configuration of the two neighboring 
sites by 

4 = ~Ji((l + m?)(l + m*)> 

+ ±(J 2 -J 3 )(l- m ?m<*) + j AF , (4) 

where mf = nf — nf are the orbital operators, and the Jaf term comes from the magnetic 
superexchange between the nearest neighboring local spins. It is worthy of pointing out here 
that the orbital operators introduced above has unusual operator algebra, quite different 



from that of the spin operators. It can be shown that they satisfy the following relations: 
(mf) 2 = l,mf + m% + mf = 0, and [m^rrv-] = [mf,m*} = [mf,raf] = \/2>(d\^d iz - d\ z d i2 ). 
The orbital Hamiltonian H Q can be written as 

Ho = Y. u tJ m?m« - £ ^ + f E ^ 
(y) (y) * 

~gJ2Qi ^ 2 cos0i + -^=(mf-mf)sin0^ , (5) 

where the effective orbital coupling depends on the spin configuration of the two neigh- 
boring sites by 

Uij = i(Ji -J 2 + J 3 )(S l ■ Sj) + (3J 3 - 2Jx + 2J 2 ) , 

and /ijj = — \ Ji(Sj • Sj —4) . All these coupling parameters J^, w^- and fojj in i7 5 and if are 
determined not only by the spin and orbital configurations of the nearest neighboring sites 
i and j, but also by the direction of the (ij) bond. For short, we denote them by J a , u a and 
h a thereafter. If there are two symmetric directions in the system, e.g., x- and ^-direction, 
one has J x = J y , u x = u y , and h x = h y . 

The spin Hamiltonian H$ is an anisotropic Heisenberg Hamiltonian with SU(2) sym- 
metry. At low temperatures, the spin configuration along the a direction is determined by 
the sign of J a . Dividing the system into two sublattices A and B according to their spin 
alignments, and performing the well-known Holsten-Primakoff (HP) transformation in the 
linear spin wave (LSW) theory, up to the quadratic terms, we diagonalize H s as 

H s = Ek(4^k + xUk + 1) - 12W] . (6) 

k 

Here ip^ and Xk are the quasiparticle operators of the spin wave excitations with k 
the wave vectors of one sublattice. The quasiparticle spectrum is given by = 
y/(W + P k ") 2 - (P^ ) 2 , with P k T = 2SJ2 a Ja cos k a e(TJ a ), and W = 2SJ2 a \Ja\, in which 
is the unit step function. 

The orbital Hamiltonian Ho looks quite like Hs, where the orbital operator may be 
regarded as an isospin operator. But the absence of the SU(2) symmetry in H Q and the 
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abnormal algebra of orbital operators make the orbital operators quite different from the 
spin operators. For example, orbital F-type arrangement is not an eigenstate of Hp, and 
in case of orbital AF configuration, on orbital sublattice A or B there are only several 
preferable orbital alignments at which the ground-state energy of the system reaches it 
minimum, unlike in a AF spin system where all the spin orientations on a sublattice are 
energy-degenerate. In this case, the orbital state at site i can be generally expressed as 
|z) = cos(9 a /2)\z) + sm(8cr/2)\z) with a = + for i G A and — for B, respectively. From the 
symmetry of h x = h y and relation + rriiy + = 0, the second term on the right-hand 
side of Eq. (5) can be rewritten in a more intuitive form 

H z = -e z £ mt , (7) 

i 

with e z = h z — h x . This anisotropic Hamiltonian arises from anisotropy of electronic hopping 
integrals in orbital space as well as unusual algebra of orbital operators. Both u a and h a 
are anisotropic and depend on the spin configurations along the a direction, as shown in 
their expressions below Eq. (5). Since J\ — J 2 + J3 is always positive Jll we have u z < u x 
and e z > for the A-type AF spin configuration; u z > u x and e z < for the C-type AF 
one; and u z = u x and e z = for the ferromagnetic (F) one. The static JT distortions 
Qi are approximately treated as classical variables and assumed to be different in the two 
sublattices, i.e., Qi = Q a and fa = cf) a with o = + (— ) for % G A (B). From x-ray diffraction 
experiments, it has been confirmed that the MnOg octahedron is elongated along the x or y 
direction, and the octahedrons are alternatively aligned in the x-y planeji! which corresponds 
to cf) + = 2n/3 and 0_ = —2it/3 in the present formula. Similar to the treatment of the spin 
degrees of freedom, we perform the HP transformation for localized orbital operators.! To 
the lowest order, H Q can be diagonalized as 

#o = E eUUk, + ~ - Pa) + E c . (8) 

Here and £k<r are the quasiparticle operators of the orbital excitations, the second term 
stands for the quantum fluctuation energy where 
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P a = - 4w a cos 61 cos Q°_ + 2e z cos 9 a 

a 

+ COS 2 (6^ - (j) a ) . 

K 

with 9% = 9 a — ip a , and Eq is the classical grand-state energy. The expression for E c depends 
on the orbital configuration. For both G- and C-type AF configurations, it is given by 

1 K 

E c/N = ^tt a cos0+ cos 6°L - -^[e^cos^ + gQ a cos(9 a - <j> ff ) - —Q 2 a ], 

with iV the number of the sites. In principle, 9 a and Q a in Eq.(§) should be determined by 
minimizing the total ground state energy of the system. In the present case, the quantum 
fluctuations in Hs and Hq are small, and so the ground-state energy can be approximately 
replaced by E c . It is found that besides the same ground-state energy E c , there is the same 
excitation spectrum for the C- and G-type AF orbital configurations, yielding 

s ka = sj\{Pl + P 2 + a[{Pl - P 2 ) 2 + 16P + P_C2]V2} , 

where C*k = E Q 2it a sin 9^_ sin cos k a . This degeneracy of C- and G-type AF orbital config- 
urations agrees to Mizokawa and Fujimori's result. @ Independent of the magnetic structure, 
such a degeneracy suggests the possibility of a mixed C- and G-type AF orbital configuration 
in the system, i.e, neighboring orbital states along the z direction may be either "parallel" 
or "antiparallel" . In the absence of the Coulomb interactions, a C-type AF orbital structure 
may have lower energy! 

The JT coupling plays an important role in determining the orbital ordering. In the 
absence of the JT coupling and in the small limit of e z , the e g electrons may occupy two 
"antiparallel" states in the two sublattices: (\z)±\z))/\/2 (9 + = — 6L = 7r/2) for u z < u x ; \z) 
and |^) (9 + = 0, 9- = —it) for u z > u x . Such symmetric "antiparallel" states will be broken 
by the uniform crystal field appeared in Eq. (7). Furthermore, the JT distortions also lead 
to an effective anisotropic crystal field acting on the two sublattices. To distinguish it from 
the uniform crystal field e z , we call it the JT field. The JT field, whose strength increases 
with the coupling constant g, tends to align the orbital states in the two sublattices towards 
\y) (9 + = 2n/3) and \x) (9- = —2tt/3), respectively. 
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The orbital ordering is described by the average value of operators mf. From the orbital 
spectrum, it can be shown that 



The second term on the right-hand side of Eq. (10) comes from the quantum and thermal 
fluctuations. To keep a good approximation, this term must be small at low temperatures. 

We now discuss the ground state of the system. First, it is impossible to realize an 
isotropic orbital ordering. Since mf + m\ + mf = 0, if (mf) = (mf) = (mf), they must be 
equal to zero and there is no any orbital ordering. From Eq. (4), it then follows that the 
anisotropy in (mf) will lead to anisotropic J a . At zero temperature, M a = 1 and (m") = 
cos#" if the quantum fluctuation in Eq.([5]) is neglected. Taking into account the symmetry 
requirement of (mf + mj) = (m\ + mj), we get two possible relations: (I) 6 + + 6_ = or (II) 
6 + — 6- = 7i. As the quantum fluctuation is taken into account, relation (I) keeps unchanged, 
while relation (II) is satisfied only approximately. In both cases, we have J x = J y ^ J z 
from Eq. (4), provided the small quantum fluctuations are neglected. Since the magnetic 
structure at zero temperature is determined by the sign of J a , the same sign of J a , regardless 
of anisotropic magnitude of them, will lead to a F or G-type AF spin configuration, while 
different signs of J x and J z will result in a A- or C-type spin configuration. Our calculations 
show that the ground-state magnetic structure is very sensitive to the on-site Coulomb 
interactions. Even though the magnetic superexchange Jaf is fixed and the JT coupling is 
absent (g = 0), an evolution of spin configuration in the order ofF—>-A—>-C^G can be 
obtained with increasing the Coulomb interactions, as shown in Fig. 1(a). It is found that 
spin configurations A and G satisfy relation (I), and spin configuration C satisfies relation 
(II). Figure 1(b) shows that an increasing JT coupling narrows gradually the C-type AF 



(m%) = M a cos6\ 



(9) 



with a = + (— ) for i G A (B), where 




(10) 
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region. This is because the JT coupling tends to align the orbital states along \x) and \y), 
and so raises the effective ferromagnetic coupling in the x-y plane and the AF coupling in 
the z direction, making the C-type AF spin configuration unstable. 

We next discuss the orbital excitation spectra. Owing to the absence of SU(2) symmetry 
in the orbital Hamiltonian, an orbital excitation spectrum usually has an energy gap. For A-, 
C- and G-type AF spin configurations, there is always an energy gap in the orbital spectrum, 
regardless whether or not the JT coupling is taken into account (not shown here). However, 
if the JT coupling is absent, gapless orbital spectra may appear for the F spin configuration, 
Furthermore, if relation (II) is satisfied, the orbital spectrum has a two-dimensional form: 
£ko- = QUxJl + cr(cos k x + cos k y )/2. For such a two-dimensional spectrum, quantum and 
thermal fluctuations, characterized by the second term of M G in Eq. (|T0|) , will completely 
destroy long-range orbital ordering at finite temperatures,^ resulting in an orbital-liquid 
state similar to that obtained by Ishihara et a/.El The orbital excitation gap can be widened 
by the JT field acting on the orbital states. It is very similar to an anisotropic magnetic 
crystal field on the spin states in an AF Heisenberg Hamiltonian. Quantum fluctuations are 
greatly suppressed by this JT field, making the orbital ordering stable. 

At finite temperatures, (Sf) and M a in Eq. (9) serve as the spin and orbital order 
parameters, respectively. Both of them decrease with increasing the temperature, and (S*) 
(M a ) vanishes as the temperature is increased beyond the critical temperature T/v (To)- 
One may evaluate and M a from a self-consistent equation for (S-) and Eq. (10). In our 
calculation, parameters Ji, J2 and J3 are taken from the Racah parameters^! and t = 0.41 
eV.0 The system is found to have an A-type AF spin configuration at low temperatures. In 
Fig. 2 we plot the variation of T/v and Tq as functions of the strength of the JT coupling. 
Both T/v and To increase with the JT coupling, but there are different physical origins. The 
increase of T/v is attributed to an enhancement of the effective magnetic coupling J x and 
J y caused by the JT field. On the other hand, the increase of To stems from the fact that 
a stronger JT field will widen the energy gap of the orbital excitation spectrum, and so a 
higher temperature is required to excite orbital quasiparticles to break the long-range orbital 
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ordering. According to experimental data and theoretical analysis, g is of the same order 
of magnitude as t and K is greater than g by a factor of ten to hundred, EE so that g 2 /K 
is the order of O.Olt ~ O.lt. According to Fig. 2, to fit with To = 780K measured by the 
experiment, the strength of JT coupling should be g 2 /K = 0.045i, at which the calculated 
T/v = 146K is very close to the experimental value of T/v = 140K. The present calculation 
may overestimate the critical temperatures due to neglecting the frequency-softened effect 
for the excitation spectrum at high temperatures, and so the required strength of JT coupling 
may be greater than the evaluated magnitude. 

This work was supported by a CRCG grant at the University of Hong Kong. 



9 



REFERENCES 
1 J. B. Goodenough, Phys. Rev. 100, 564 (1955). 
2 Y. Murakami et aL, Phys. Rev. Lett. 81, 582 (1998). 
3 K. I. Kugel and D. I. Khomskii, ZhETE 64, 1429 (1973). 
4 S. Ishihara, J. Inoue and S. Maekawa, Physica C 263, 130 (1996). 
5 A. Moreo, S. Yunoki and E. Dagotto, Science 283, 2034 (1999). 
6 R. Shiina, T. Nishitani and H. Shiba, J. Phys. Soc. Jpn. 66, 3159 (1997). 
7 R. Maezono, S. Ishihara, and N. Nagaosa, Phys. Rev. B 58, 11583 (1998). 

8 J. van den Brink, P. Horsch, F. Mack and A. M. Oles, Phys. Rev. B 59 6795 (1999). 

9 T. Hotta, S. Yunoki, M. Mayr and E. Dagotto, |cond-mat/9907034| (1999). 

10 A. J. Millis, Phys. Rev. B 55, 6405 (1996). 

11 S. Q. Shen and Z. D. Wang, Phys. Rev. B 59,3291 (1999); |cond-mat/9906l26| (1999). 

12 L. F. Feiner and A. M. Oles, Phys. Rev. B 59, 3295 (1999). 

13 A projective perturbation approach yields J x = t 2 /[8(U + 3J H /2)}, J 2 = t 2 /[10(U' + 3J/2 + 
3Jff/2)], and J 3 = t 2 /[10(U' - J/2)]. Here t is the hopping integral, U (W) is the intra- 
(inter-) orbital Coulomb interaction, J and Jh are the Hund's coupling between the e g 
electrons and between the e g and the t 2g electrons, respectively. 

14 J. Kanamori, J. Appl. Phys. 31, 14S (1960). 

15 G. Matsumoto , J. Phys. Soc. Jpn. 29, 606 (1970). 

16 T. Mizokawa and A. Fujimori, Phys. Rev. B 51, 12880 (1995); ibid 54, 5368 (1996). 
17 P. C. Hohenberg, Phys. Rev. 158, 153 (1967). 

18 S. Ishihara, M. Yamanaka and N. Nagaosa, Phys. Rev. B 56, 686 (1997). 

10 



U, U' and J are connected through three Racah parameters by U = A + 4B + 3C, [/' = 
A - B + C and J = 55/2 + C. Taking A w 4.7 e y [Ref. 12], B = 0.107 and C = 0A77eV 
[A. E. Bocquet et al, Phys. Rev. B 46, 3771 (1992)], one has U ss 6.6, C/' w 5.1 and 
J = 0.745eV. 1.2eV is estimated from Lax.^Sr^MnOa at x=0.175. [M. Imada et al, 

Rev. of Mod. Phys., 70, 1039 (1998)]. 

Y. Okimoto et al, Phys. Rev. Lett. 75, 109 (1995); Phys. Rev. B 55, 4206 (1997). 

A. J. Millis, Phys. Rev. B 53, 8434 (1996); Y. Motome and M. Imada, J. Phys. Soc. Jpn. 
68, 16 (1999). 



11 



FIGURES 

FIG. 1. Phase diagram at zero temperature in the absence (a) and presence (b) of the JT field. 
The parameters used are Jaf = 0.006 and Jh = 4/3 with t the unit of energy. The relation 
U = V + 2 J has been usedlH and U = 20 fixed in (b). 

FIG. 2. Critical temperatures Tjv and Tq as functions of g 2 / 1 Kt. 
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